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Быстрый алгоритм решения задачи 
о назначениях для нахождения нижней 
границы стоимости маршрута коммивояжера 


Показано, что если в оптимальном решении задачи о назначениях (ЗН) и ее матрице стоимостей порядка 
п заменить значение какого-либо элемента на бесконечно большое число, то оптимальное решение ЗН 


ДлЯ полученной матрицы находится за время 0 (”) й 


Цель работы и постановка задачи 


Сформулируем задачу о назначениях (ЗН) для матрицы стоимостей С = Е | ‚ где 
с. =® при =] и с,е2у или с, =® при 1*], Ь]=Ьп, 2 — множество 


неотрицательных целых чисел. Найти 
п 
С( с) = тт У` с, 
= 
здесь л= (=[|] : п[2] о л["]) — перестановка и столбцов матрицы, а 
а = (< , с|[2],...9[п]) — оптимальная перестановка, или решение ЗН стоимостью 


п } РИ > 

С( о) р = ме 5 бой #00, 1=1, и. Заметим, что для матрицы С, содержащей 

элементы с, =, 1 /, ЗН может не иметь решения. Предположим, что получена 
перестановка ос’. 

Рассматриваемая здесь задача состоит: а) в построении по известному решению 

с решения о, ЗН с матрицей стоимостей С, , отличающейся от С только тем, что 


„; 6) время нахождения о, меньше, 


элемент с, е [бе с.) в С равен © в С 
чем время, за которое можно определить о’. 
Матрица С совпадает с матрицей стоимостей М№Р-трудной гамильтоновой 
задачи коммивояжера, решением которой является циклическая перестановка 
минимальной стоимости. Поскольку ЗН эффективно разрешима, то она используется 
как релаксация в точных и приближенных алгоритмах решения №Р-трудных задач 
проблемы коммивояжера. Точные переборные методы построения циклических 
перестановок характеризуются обращением к уже полученному решению ЗН о’ для 


нахождения нового решения в матрице, где один из элементов Сой заменен на ©, 


ТЕ {1, и п}. К ним относится алгоритм Кристофидеса, который строит по схеме 
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метода ветвей и границ маршрут коммивояжера в полном взвешенном графе, 
заданном симметричной матрицей стоимостей [1]. Алгоритм Кристофидеса начинает 
работу с решения ЗН, позволяющего оценить снизу стоимости всех возможных 
замкнутых маршрутов и переходит к выбору в полученном решении элемента, 
инициирующего ветвление. После замены значения выбранного элемента на с 
снова находится решение ЗН, по которому определяется нижняя граница стоимости 
маршрутов коммивояжера, не включающая этот элемент. 

'Трудоемкости всех известных алгоритмов решения ЗН оценивается величинами, 


не меньшими 0(" [1], [2]. Поэтому неоднократное обращение к какому-либо из 


них, характерное для метода ветвей и границ, требует значительных временных затрат. 
Эти затраты могут быть уменьшены нахождением перестановки о,, алгоритмом, 


временная сложность которого меньше, чем 0(»") . В его основу положены подходы 


к решению задач о паросочетаниях [1], [3]. 


Сведение к задаче построения кратчайшего пути 


Матрице стоимостей ЗН С поставим во взаимно-однозначное соответствие 
двудольный граф (х эр ) : 


Хх | — |у | =и, где вершина 1е Х соединена с вершиной 
ТЕУ ребром [1, ЕЦ весом с; *®. 


Паросочетанием графа называется такое множество ребер, что никакие два из них 
не имеют общей вершины. Ребра графа, не входящие в паросочетание, называются 
свободными. Максимальное паросочетание — паросочетание с наибольшим числом 
ребер для данного графа. Вершина, принадлежащая ребру паросочетания, называется 
насыщенной, остальные вершины графа — свободными. Если Г — насыщенная вершина 


ребра [2, Л ‚ входящая в паросочетание, то вершина / определяется, как напарник 1. 


Совершенное паросочетание — паросочетание, насыщающее все вершины графа. 
В двудольном графе (Х,У,И), соответствующему матрице С, мощность совер- 


шенного паросочетания, если оно существует, равна п [2]. 
Решение ЗН о является совершенным паросочетанием двудольного графа 


(Х,У,0), включающим ребра с минимальным суммарным весом. Элемент (х, у) 
матрицы С представлен в паросочетании с’ ребром [х, у], которому присваивается в 
матрице С, вес, равный со. Изменение значения с„, © 70 на о влечет удаление из 
с ребра [х, у]. Таким образом, поставленная задача заключается в нахождении 
совершенного паросочетания о, с минимальным суммарным весом С (о) 
относительно заданного паросочетания ©—|[х, у]. 

Рассмотрим чередующийся путь (у,ъ,,....у,...у„) графа (Х, УИ -[х, У] , 
УЕХУУ, 1=Ьт, т. е. путь, в котором ребра (У, 2), (Уз Уд}... (Уже-ь Ук } 2. 
свободны, а ребра [у,%], [4,%],.... [Ужин |»... Входят в паросочетание. 


Чередующийся путь и, у›,...› Уи | называется увеличивающим, если вершины у и 
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у„ свободны [2]. Увеличивающий путь состоит из внешних и внутренних вершин, 


расположенных соответственно на нечетных и четных позициях. Поэтому в нем 
первое и последнее ребра свободны. 


Лемма 1. Решение о., существует, если в графе (х ‚УИ -[х, У], содержащем 
паросочетание о-[х, у], существует увеличивающий путь из вершины хеХ в 
вершину уЕУ. 

Доказательство. Заданное паросочетание о-[х,у| состоит из и-1 ребер. 
Следовательно, в графе (х ‚У, (- [8 У] только две вершины х и у свободны. Если в нем 
существует увеличивающий путь Р между вершинами х и у относительно паросочетания 
о-[х, у], то существует паросочетание (с —[х, У] ФР= ((с =[х, У] -Р) о(Р- 
-(< —[х, у] мощностью и, а значит, и совершенное паросочетание о’... 

Лемма 2. Решение ЗН о, состоит в нахождении кратчайшего увеличи- 
вающего пути Р., из вершины х в вершину у относительно паросочетания ©’ [х, у]. 

Доказательство. Пусть в графе (Х,У,И -[х,у]) построен кратчайший увели- 
чивающий путь Р., из вершины х в вершину у относительно паросочетания ©’—|[х, у]. 


Определим совершенное паросочетание о’, = (( -[»>])-2. ) © (Р -(<-[х У] 


х 
и покажем, что оно имеет минимальный суммарный вес ребер С (о) среди всех 
совершенных паросочетаний графа (Хх эй ОГ [х, У]. 
Рассмотрим о, вместе с совершенным паросочетанием о’, построенным в 
графе (х ‚У, Ы ) . Очевидно, для о’ и о,, справедливо неравенство С (с) 8 (<„) : 


Числа С (<) иС (о) содержать общую часть 


С(ху)= > с. 


<[Йеохуп о 

Из неравенства С(о’)< С( о.) для 9 По = (с — [х, У] = | 2:2 )- [х, У] 

следует неравенство 
С(о)-С(х, У) = С(о)-С(х, У). 

Но С (о„)-С (% У) — минимальная сумма весов ребер множества Р,, - (< - 
—[х, У] ‚ поскольку Р,, — кратчайший увеличивающий путь из вершины х в вершину 
у относительно паросочетания ©’— [х, У] ) 

Лемма 3. Если в графе (Хх „У; ) содержится и совершенных паросочетаний 
°„,и>2, включающих ребро [х, у] и имеющих минимальную стоимость С(о’ ), 
т = ПА, ‚ то длина кратчайшего увеличивающего пути в графе (Хх ГЫ = [х, У] не 


зависит от того, относительно какого из паросочетаний бт [х, У] построен этот путь. 
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Доказательство. Без потери общности предположим, что граф (Х,У,И) 


содержит два оптимальных совершенных паросочетания 0\ = (с [1], ва 


ол [п] и в, = (95 [|], 05 [2],...9>[п]), С(в,)=С(9>), [х, у] 1 оо. 


Так как о\ * 0), то 0% = (о: -[х “(> -[ у], (о-в) “(о-в )=9, 


|1 = бу = |б> —о0|, С(о' _ бо) = С (с. —: бо) . По отношению к подмножеству ребер 
о! - о паросочетания о, -[х, у| каждое ребро о.[е о. -—0% является свободным, 
и обратно, каждое ребро о[Пео1-©у свободно относительно подмножества 


о) —оу паросочетания д. =[х, У] (рис. 1). Кратчайший увеличивающий путь Р.) из 


у 
вершины х в вершину у относительно паросочетания ©-— [ь У] включает 
подмножество ребер о1—0% и подмножество свободных ребер о.-—о%у и имеет 
длину 2С (о = бу) . Но такую же длину имеет кратчайший увеличивающий путь Р,,› 
относительно паросочетания о’, - [х, у], содержащий подмножество ребер ©, —0% и 
подмножество свободных ребер о — [х, у]. 

Следствие. Пусть о„ — совершенное паросочетание минимальной стоимости в 
графе (х из ) ‚т= ПТ, ‚ Ш22. Стоимость совершенного паросочетания с’„, графа 


(Хх ‚У, Ц = [х, У] не зависит от паросочетания о„ =[х, У|, относительно которого 
построен кратчайший увеличивающий путь из вершины х в вершину у в графе 
(х,у,и-[х,у]}. 

Доказательство. Достаточно положить д=2. В результате построения 


кратчайших увеличивающих путей Р. 


я и Б.,> из вершины х в вершину у в графе 


(х, УИ -[х, У] относительно паросочетаний 01 -—[х,у| и о. -[х, у] получим два 
совершенных паросочетания о.) И ©.у2, каждое из которых включает подмножество 
ребер оу = (в -[х, у]^о> —[х, у]) . Определим их стоимости: СС) =С(°%)+С(в.- 


—0), С(о) = С(°,)+С(01 0%) ‚ но С(о 0%) = С(с. 0%). 
Доказанные утверждения дают основание перейти к организации эффективного 


поиска О с помощью нахождения тех увеличивающих путей из вершины х в 


вершину у относительно паросочетания о’, =[х, У], среди которых содержится 
кратчайший путь Р,,. Поэтому для решения поставленной задачи применим адапти- 


рованный к ее условиям метод поиска в ширину, развивающий идеи алгоритма 
построения максимального паросочетания в двудольном графе [2]. 
Как и в алгоритме из [2], отнесем х к внешним вершинам искомого пути Р,. 


Тогда у — внутренняя вершина в . В двудольном графе (х оу = р У] ‚ содержащем 
паросочетание о’ — [№ У] ‚ все вершины Те Х и /Е У за исключением х и у насыщены. 


Будем искать путь Р» построением всех возможных чередующихся путей из 
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вершины х. В каждой паре ребер ((&., Я ( У , чередующегося пути с началом в 
вершине х вершины й,,1„еХ внешние, а вершина д еЕУ - внутренняя. Так как в 
нем вершина х свободна, то свободно ребро (х, Л), а ребро (Л,#„) содержится в 


паросочетании ©’, =[х, У]. 


у У. У Уи = У Ур =, № у. № т 


У хгУ, У; № о О а 9: И 


Рисунок | - 0% = [м2], [м] а) 1—0 [5 4] [у 6 | [8] 


6) ©› —0% = [5%], [№74], [№9 %]} 


Идея эффективного поиска пути Р„ в ширину заключается в пошаговом 
построении орграфа, в котором вершина х имеет нулевую полустепень захода, а 


вершина у — нулевую полустепень исхода [3]. Каждая пара ребер ((х, Л), (ль и), 


представлена в орграфе дугой [8 в) весом с =с,, +с,„. Поскольку й, является 


1 ЕЛ 
напарником /, а любая вершина не может иметь более одного напарника, то такая 


замена однозначно определяет очередную внешнюю вершину чередующегося пути, 
игнорируя внутреннюю вершину. Вершины й, ‚ инцидентные х, образуют первый ярус 


вершин орграфа. Тогда для каждой вершины й,, смежной со свободной вершиной у, 


можно построить увеличивающий путь (х, Л» й У), определить соответствующее 


ему совершенное паросочетание [х, И, У (в -[х У|-[л]), вес которого 
п 

равен с) +с,,+ + —‹„-с,). Среди всех полученных совершенных 
ри 


Е 1 1 
паросочетаний выберем паросочетание о’, с наименьшим весом С (®»,) ‚ принимаемым 


за рекордное текущее значение МИМ верхней границы для С (25) . В каждом 
построенном увеличивающем пути отметим дугу, входящую в у, и вершину 1, , если 


она не смежна в (х ИО =[х, У] с внутренними вершинами или смежна только с 


вершиной у. Паросочетание в. 


з тогда является решением поставленной задачи, 


когда оно имеет вес, равный М/Х, и все вершины й, ‚ а следовательно, и дуги (1, У) 


отличны (рис. 2). 
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Рисунок 2 — а) Двудольный граф (х Е [х, ы] и его паросочетание 
а-[х, У] = {[1, 2] [6,3], [5, м} [4 В 6) увеличивающие пути из вершины 
х=& ввершину у =й относительно паросочетания о’— Е: У] ; в) двудольный граф 


(х, УИ -[х, У]) и его совершенное паросочетание о’, = {[, ВЫ л] 
[54,5] 


В общем случае каждый ярус вершин орграфа, за исключением непосредственно 
предшествующего у, состоит как из отмеченных, так и неотмеченных вершин. Чтобы 


найти вес вершины второго яруса, для каждой неотмеченной вершины й&, определим 
ее вес С(й)=с, [Ж‚ Леа -[х, у]. Не теряя общности, предположим, что в графе 
(Х, УИ -[х, У] вершина й, инцидентна двум свободным ребрам (&,, /.) и (1, /.}, 
вершина {, — одному свободному ребру (1„, /,), 1 *у, Л #у, весом С(,)= су, 


[2, т] Е а-|[х, У] ‚ г, является напарником #., а 1, — напарником /,. Представим пары 
ребер ((&,1,),(#51,)), (и, /.), (ль), (ыы л)(»-и)) графа (Ж,У,И -[ж]) 


соответствующими дугами (^., 1 ал (+, т 5 в) в строящемся орграфе и отнесем 
Г, › & ко второму уровню его вершин (рис. 3). 


Вычислим веса И И и: 


С(е,) = С(и) че, 1, [Ев -[х, |; 
С(ь) =тт Ка) +» СЕ.) Сил } у |, М Е а-|[х, У] ь 
Пусть 1, — вершина второго яруса, из которой исходит дуга (5 У) . Вершина (1, 
в графе (хХ „зе [8 У] принадлежит увеличивающему пути из х в у, длина которого 
равна С ее Предположим, что это путь (х, Л, и, №,» У). По нему найдем со- 


вершенное паросочетание о’, = (> Ям» (9 -[ь №14] 
с суммарным весом ребер 


п 
С(о р ыы Са» ео Ус м о Е 
15 


«Штучний 1нтелект» 42011 411 


Левченко А.Ю., Морозов А.В., Панишев А.В. 


отметим дугу (%.>) и вершину Г,, если она не смежна в (х, УИ -[х, |) ни с 
какой другой внутренней вершиной. После выполнения таких действий для каждой 


вершины второго яруса, инцидентной у, определим совершенное паросочетание сб, 


с минимальным суммарным весом ребер С (25, и, если С (©. < МИ ‚ то положим 


му= = С(о>,}. 


т У 5 | 
Рисунок 3 — Увеличивающие пути из вершины х, 


содержащие дуги (х,1,), (х, и), (№), (т а 


Процесс нахождения увеличивающих путей продолжается для неотмеченных 
вершин яруса, сформированного на очередном шаге, и завершается построением 
орграфа, в котором отмечены все дуги, входящие в вершину у. Искомым решением 


с. является полученное на шаге и, и=1, 2,..., совершенное паросочетание О. с 


суммой весов ребер, равной М/МУ. 


Поиск ©, в ширину упрощается в орграфе (У, А), где И={м,у,,... у}, 
| = 74 = |у ‚м, =Х, у, =у, А - множество дуг, таких, что (увьу,) = А, если и только 


если в двудольном графе (х, УИ -[х, У]: а) (&, Л) И-[х, у] и л=У; 6) вершина 
1ЕХ смежна с напарником ]„еУ вершины 1ЕХ, т.е. вершины й,]„,й 
образуют пару ребер (&,/„), [/„›й| некоторого увеличивающего пути. В орграфе 


с, ‚ иначе. Нетрудно видеть, 


(У, А) дуга имеет вес с и бы 


им =Сил› ®6ли Д=у, ис 
что решение о’„, достигается построением во взвешенном орграфе (У, А) кратчайшего 


пути из вершины у. =х в вершину у. =у. 
т уу, -У 


Алгоритм решения задачи 


Изложенную схему эффективного поиска ©., представим алгоритмом, выполня- 


ющим следующие действия. 
50. Алгоритм нахождения решения ЗН о., для матрицы стоимостей С, по 


известному решению ЗН о для матрицы стоимостей С; С — матрица порядка п, 


элемент которой с; =® при {=] и се или с, = при {*], Б]=Ьп, 2 = 
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множество целых неотрицательных чисел; Су отличается от С тем, что в С 
бе Е: Са су=°®, ха у; су, является одним из элементов ст # о, Ем, В 
решении © = = (< 1), о|[2], уе <["]) , 

$1. По матрице С„, и паросочетанию а-[х, У] построить взвешенный орграф 
с множеством вершин Г = У, аа у, } и множеством дуг А следующим образом: 


для каждого элемента (К, 7) матрицы С 


у› Су #00, соответствующее ребро которого 


не входит в паросочетание ©’-[х, у|, если /*у, то (у, у) ЕЛ, с в 


у тогда и 


У > 
только тогда, когда элемент (1,/) представлен ребром [1,7] в нпаросочетании 


а-[х, У] ‚ иначе (%ь у, ЕЛ, Ст; = Су. 

52. В построенном орграфе (т, А) алгоритмом Дейкстры найти все кратчайшие 
пути из вершины у,. 

53. Если в орграфе (7 А) не существует пути из вершины у, в вершину у, , то 
конец: ЗН для матрицы Су не имеет решения. 

534. Каждую дугу (Ук, у.) кратчайшего пути из вершины у, в вершину у, заменить 

. . ! 

на ребро [К, /| весом су; сформировать из полученных ребер подмножество 0”, 
искомого паросочетания о... 

$5. Найти все ребра паросочетания о.,, выполнив такие действия: исключить 

1 

из паросочетания о-[х,у| каждое ребро [А,г"], если о’, содержит ребро [(, $], 


г=5; объединить подмножество из оставшихся ребер множества о—[х, у] с бу; 


вычислить С [е5) 3 
Теорема. Представленный алгоритм корректно находит за время 0(») решение 


ЗН 9’, для матрицы стоимостей порядка п. 


Доказательство. Алгоритм останавливается на шаге $33, если в построенном 
орграфе (т, А) не существует пути из вершины у, в вершину у, . В этом случае по 


лемме | не существует увеличивающего пути из вершины х в вершину у в 
двудольном графе (х ‚УИ -[х, У] ‚ соответствующем орграфу (У, 4). Следовательно, 
граф (Хх орз О —[х, У] не содержит совершенного паросочетания, и поставленная 


задача не имеет решения. 
Действия по построению орграфа (У, 4), выполненные на шаге $54, требуют 


времени, оцениваемого величиной О(»?). Такую же временную сложность имеет 


шаг $52, на котором алгоритмом Дейкстры находится кратчайший путь из вершины у, в 


вершину у, . Шаги 54 и 55 имеют одинаковую оценку трудоемкости, равную О(п) } 
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Пример. Для матрицы 


с 35 < 14 < 5| ®х 18 
со © 14 ® 120 41 51 93 


64 [2] 51 ® 97 ® ® ® 
[61] ю ® ® 22 ® ® ® 


64 ® 8 ююоохоою 


перестановка @= [7 6, 8, 5, 4, 2,1, 3) является решением ЗН; С (о) = 7+ + 


о мин офЬь- 
8 
> 
> 
8 
в 
= 
8 
8 
8 
8 


+ С3з + С4з + 654 + сво + си + сз = 61+12+18+20+20+12+61+18=222. Определим матрицу 
С.5 ‚ заменив в С значение с,; =20 на с. Выполним алгоритм нахождения решения 
045 для матрицы Сдз. 

Матрице С.5 соответствует двудольный взвешенный граф (х ‚У, Ц - [4, 5] ‚В 


котором найдем кратчайший увеличивающий путь из вершины х=4 в вершину у=5 
относительно паросочетания а- [4,5] (рис. 4). 


На рис. 5 изображен граф (У, А) ‚ построенный на шаге 31 алгоритма нахождения 
045 . В орграфе (У, А) из вершины 1 исходит дуга в вершину 8 ‚ поскольку с; ©, а 
элемент (8, 3) матрицы С.5 соответствует ребру [8, 3] весом 18, входящему в 
паросочетание о-|4, 5]. Так как у=5 и с15 *°, то орграф (У, А) содержит дугу (1, 


5) весом, равным с15 =59. 


Рисунок 4 — Граф (Х, У, -[4, 5] ; 


утолщенными линиями обозначены ребра паросочетания ©’ — [4, 5] 
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В результате выполнения алгоритма Дейкстры определим кратчайший путь в 
каждую вершину графа (У, 4), достижимую из вершины 4: 


1 р 3 4 5 6 7 8 
51 41 92 0 96 128 78 14 


4(93) — 6(51) 3(47) 1(64) 1(64) 97 


Рисунок 5 — Взвешенный орграф (У, А); утолщенными линиями обозначены дуги 


кратчайшего пути из вершины 4 в вершину 5 


Кратчайший путь из вершины 1 в вершину 5 в графе (У, А) состоит из дуг (4, 8), 
(8, 7), (7, 5) с весами соответственно сз =14, с =64, с; =22 (рис. 5). В графе 
(Хх м - [4,5] ему отвечает кратчайший увеличивающий путь Р.5 = {(ы, 23), (1в), 
(в, л).(л.5).(5,/5)}, в котором ребра (и, Л), (&,Л), (&,/5) являются сво- 
бодными. Они образуют подмножество ребер 0/5 = {[4, 3], [8,1], [7,5] искомого 
паросочетания 0.5. Чтобы определить все ребра в 0.5, нужно в паросочетании 
о-[4,5] = [17], [7,1], [2, 6], [6,2], [3,8], [8, 3], [5,4]} исключить ребра [8,3] и [7,1] 
и объединить полученное подмножество с подмножеством ребер {[4, 3], [8,1], [7, 5] | 
В результате искомым решением является паросочетание од; ЕТ [7,5], [5,4], 


[4,3], [3,3] [8,1], [2,6], [6,2]} стоимостью С(945) = 233. 


Выводы 


Изложенные результаты получены в процессе обсуждения и поиска возможно- 
стей, уменьшающих трудоемкость известных алгоритмов типа ветвей и границ для 
нахождения кратчайших маршрутов коммивояжера. Как известно, эффективность ме- 
тода ветвей и границ зависит от способа получения оценок, ограничивающих снизу в 
выбранной схеме ветвления стоимость оптимального решения комбинаторной задачи. 
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Для задач класса коммивояжера существует несколько способов нахождения оценок. 
В данной статье утверждается, что способ с использованием в качестве релаксации ЗН 
позволяет ускорить вычисление нижней оценки стоимости циклических маршрутов со 
значением, достаточно близким к стоимости оптимального маршрута. 
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А.Ю. Левченко, А.В. Морозов, А.В. Панйиев 

Швидкий алгоритм розвязання задач! про призначення для знаходження 

нижньо: меж! вартост! маршруту комвояжера 

Показано, що якщо в оптимальному розв’язку задач! про призначення (ЗП) та й матрищ вартостей порядку 
п замнити значення якого-небудь елемента на неск!нченно велике число, то оптимальний розв’язок 
ЗП для отримано! матриц! знаходиться за час () (*) : 


А.Уи. Геусйепко, А.Т. Моготоу, А.Г. Раш5йеу 

Аоог т г Еа${ Зошйоп оЁ фе РгоШет о{ АПосайоп Гог Гоуег Глий УУео ВЕ оЕ За]езтап ВКоще 

Г 15 зпо\мп фа 1 0 тер!асе Фе уа[ае оЁ апу е!етепё Бу шйпцейу 1агое питбег ш фе орйта]| зооп Юг 
{фе рго Мет оё аЙосайоп$ ап4 Гог Из тайлх оЁ уе ИВ 517е и, феп Фе орйта! зо[авоп оЁ Фе ргоМет 


огаПосайоп Гог гези те таблх 15 опа ша ите О( п?) ь 


Статья поступила в редакцию 22.06.2011. 
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